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Abstract— We  describe  degenerate  four-wave  mixing  experi¬ 
ments  on  ZnSe  and  CdTe  semiconductor  samples  with  pico¬ 
second  laser  pulses  at  wavelengths  below  the  bandgap.  Nonlin¬ 
earities  of  third,  fifth,  and  seventh  order  are  observed  and  the 
mechanisms  for  each  are  identified.  In  all  of  our  measure¬ 
ments,  we  observe  a  fast  third  order  nonlinearity.  For  two-pho- 
ton  absorbers,  this  is  attributed  to  contributions  from  both  the 
real  (refractive)  and  imaginary  (absorptive)  parts  of  the  third- 
order  susceptibility.  Below  the  two-photon  absorption  edge,  the 
nonlinearity  is  purely  refractive.  The  higher  order  effects  are 
due  to  carriers  generated  by  multiphoton  excitation.  In  ZnSe 
at  0.532  pm,  carriers  are  generated  by  two-photon  absorption 
such  that  a  fifth  order  nonlinearity  arises  from  the  change  in 
index  due  to  these  carriers,  a  sequential  xl31:x(1>  nonlinearity. 
From  such  measurements  we  determine  the  refractive  index 
change  per  photoexcited  carrier  pair  and  the  density  depen¬ 
dence  of  the  carrier  diffusion  coefficient.  Analogous  signals  are 
observed  in  CdTe  at  1.064  pm.  The  seventh  order  nonlinearity 
observed  in  ZnSe  at  1.064  pm  results  from  the  refractive  index 
contribution  of  carriers  generated  by  three-photon  absorption. 


I.  Introduction 

WE  report  a  series  of  picosecond  degenerate  four- 
wave  mixing  (DFWM)  studies  conducted  in  ZnSe 
and  CdTe  at  wavelengths  of  0.532  and  1.064  /im.  The 
DFWM  signal  shows  a  fast  third  order  nonlinearity,  as 
well  as  higher  order  slowly  decaying  nonlinearities  due  to 
multiphoton  absorption  generated  carriers.  We  attribute 
this  signal  to  the  combined  effects  of  the  real  and  imagi¬ 
nary  parts  of  the  third  order  susceptibility  xl3)-  The  imag¬ 
inary  part  corresponds  to  two-photon  absorption  (2PA), 
while  the  real  part  is  due  to  bound-electronic  nonlinear 
refraction  (index  n2 ),  as  opposed  to  a  free-carrier  effect 
r  1  ] •  From  our  measurements,  we  obtain  the  absolute  value 
of  the  third  order  nonlinear  susceptibility  for  both  ZnSe 
and  CdTe.  This,  combined  with  independent  2PA  mea- 
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surements.  allows  us  to  extract  the  real  part  of  the  sus¬ 
ceptibility  which  corresponds  to  the  nonlinear  refraction. 

We  also  observe  a  rapid  third  order  nonlinear  effect 
which  is  dominant  at  low  incident  irradiances  in  both 
semiconductors  when  all  three  beams  are  temporally  co¬ 
incident  (zero  delay).  Fifth  and  seventh  order  nonlinear 
effects  are  evident,  depending  on  wavelength,  when  the 
gratings  are  probed  at  time  delays  greater  than  the  pulse- 
width  to  eliminate  the  signal  from  the  fast  third  order  ef¬ 
fect.  These  higher  order  refractive  nonlinearities  are  at¬ 
tributed  to  the  refractive  effect  of  carriers  generated  by 
2PA  (fifth  order)  and  3PA  (seventh  order).  Most  of  the 
experiments  discussed  here  consist  of  the  generation  of  a 
modulated  carrier  density  or  carrier  grating  created  by  in¬ 
terfering  two  of  the  three  beams  made  coincident  in  the 
semiconductor  sample.  Diffraction  of  a  third  beam  from 
this  grating  produces  a  DFWM  signal  yielding  informa¬ 
tion  on  the  nonlinearities  resulting  from  the  photogener¬ 
ated  carriers  and  their  decay  mechanisms  [2],  [3].  Decay 
of  this  signal,  which  is  the  phase  conjugate  of  one  of  the 
two  interfering  beams,  takes  place  due  to  carrier  diffusion 
and  recombination.  An  expression  for  the  diffraction  ef¬ 
ficiency  at  long  temporal  delays  is  obtained  from  coupled- 
wave  theory  [4],  This  expression  leads  to  a  calculation  of 
the  index  of  refraction  change  per  carrier  pair  per  unit 
volume  generated  via  2PA  in  ZnSe  at  0.532  pirn.  Mea¬ 
surements  of  the  grating  decay  for  several  pump-probe 
angles  gives  values  for  the  carrier  diffusion  and  recom¬ 
bination  lifetimes  in  ZnSe  at  0.532  pirn. 

After  describing  the  experimental  techniques  in  Section 
II,  we  identify  in  Section  III-A  the  dominant  nonlineari¬ 
ties  in  the  two  materials  ZnSe  and  CdTe.  We  determine 
these  to  be  fast  third  order  nonlinearities,  due  to  the  same 
processes  which  give  rise  to  the  effects  of  bound-elec¬ 
tronic  refraction  and  two-photon  absorption,  while  higher 
order  effects  are  due  to  free-carrier  refraction.  In  Section 
III-B  we  describe  our  measurements  of  the  absolute  mag¬ 
nitude  of  the  combined  third  order  susceptibilities.  Stud¬ 
ies  of  higher  order  effects  due  to  free-carrier  gratings  are 
discussed  in  Section  III-C.  In  order  to  obtain  a  quantita¬ 
tive  measurement  of  the  carrier  induced  nonlinearities,  in 
Section  III-D  we  develop  an  expression  for  the  diffraction 
efficiency  of  these  carrier  gratings  and  hence  find  a  value 
for  the  free-carrier  refractive  index  coefficient  in  ZnSe. 
By  measuring  the  angular  dependence  of  the  grating  de¬ 
cay,  we  determine  the  carrier  diffusion  coefficient  as  a 
function  of  carrier  density. 


00 18-9 1 97/91  $01. 00  ©  1991  IEEE 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

OCT  1991 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1991  to  00-00-1991 

4.  TITLE  AND  SUBTITLE 

Degenerate  Four-Wave  Mixing  Measurements  of  High  Order 
Nonlinearities  in  Semiconductors 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Central  Florida, CREOL,  Department  of 

Physics, Orlando, FL, 32816 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

We  describe  degenerate  four-wave  mixing  experiments  on  ZnSe  and  CdTe  semiconductor  samples  with 
picosecond  laser  pulses  at  wavelengths  below  the  bandgap.  Nonlinearities  of  third,  fifth,  and  seventh  order 
are  observed  and  the  mechanisms  for  each  are  identified.  In  all  of  our  measurements  we  observe  a  fast 
third  order  nonlinearity.  For  two-photon  absorbers,  this  is  attributed  to  contributions  from  both  the  real 
(refractive)  and  imaginary  (absorptive)  parts  of  the  thirdorder  susceptibility.  Below  the  two-photon 
absorption  edge,  the  nonlinearity  is  purely  refractive.  The  higher  order  effects  are  due  to  carriers 
generated  by  multiphoton  excitation.  In  ZnSe  at  0.532  f!m,  carriers  are  generated  by  two-photon 
absorption  such  that  a  fifth  order  nonlinearity  arises  from  the  change  in  index  due  to  these  carriers,  a 
sequential  x(J):  x(l)  nonlinearity.  From  such  measurements  we  determine  the  refractive  index  change  per 
photoexcited  carrier  pair  and  the  density  dependence  of  the  carrier  diffusion  coefficient.  Analogous  signals 
are  observed  in  CdTe  at  1.064  f!m.  The  seventh  order  nonlinearity  observed  in  ZnSe  at  1.064  f!m  results 
from  the  refractive  index  contribution  of  carriers  generated  by  three-photon  absorption. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


CANTO-SAID  et  a!.:  HIGH  ORDER  NONLINEARITIES 


2275 


;Nd:YAG  Laser  \ 


u  Ar 


BS2  BS3 

;di 


,:^M3 


Fig.  1.  Schematic  of  experimental  DFWM  apparatus.  0;  is  the  input 
pulse  energy  monitor,  while  D2  monitors  the  phase-conjugate  signal  pulse 
energy. 


II.  Experimental  Technique 

In  our  experiments  the  “backward”  DFWM  geometry 
is  used.  A  schematic  of  the  experimental  geometry  using 
single  43  ps  (FWHM)  1 .064  /tm  pulses,  or  30  ps  (FWHM) 
0.532  pm  pulses  is  shown  in  Fig.  1.  The  picosecond, 
Gaussian  spatial  profile  pulses  are  derived  from  a  Q- 
switched  mode-locked  Nd :  YAG  laser  system  operating 
at  1 .064  /tm.  An  electrooptic  switch  between  the  oscilla¬ 
tor  and  amplifier  ensures  single  pulse  performance.  A  sec¬ 
ond-harmonic  crystal  (KDP)  produces  the  0.532  /cm 
pulses.  This  single  pulse  is  divided  into  three  pulses 
which,  after  passing  through  variable  time  delays,  are  in¬ 
cident  on  the  semiconductor  samples.  The  three  pulses 
can  be  independently  adjusted  in  amplitude  and  polariza¬ 
tion  using  half-wave  plate  and  polarizer  combinations. 
Two  strong  beams,  forward  (Ef)  and  backward  (Eh) 
pumps,  of  approximately  equal  irradiance  are  incident  on 
the  semiconductor  from  counterpropagating  directions.  A 
weaker  beam,  the  probe  (Ep),  is  incident  on  the  sample  at 
an  angle  6  with  respect  to  Ef.  The  grating  spacing  deter¬ 
mined  from  the  angle  8  can  be  varied  from  1.2  to  8  /tm 
for  the  experiments  at  0.532  /tm.  At  1.064  /tm  the  grating 
spacing  is  fixed  at  8  /tm.  The  conjugate  wave  Ec,  which 
retraces  the  path  of  Ep,  is  detected  by  a  large  area  inte¬ 
grating  photodiode  as  are  various  reference  beams.  These 
detectors  are  calibrated  against  pyroelectric  energy  mon¬ 
itors.  All  pulsewidths  quoted  are  measured  by  autocorre¬ 
lation  using  a  second-harmonic  generator,  while  all  quoted 
spot  sizes  were  measured  in  both  horizontal  and  vertical 
directions  at  the  sample  position  by  the  method  of  scan¬ 
ning  pinholes. 

The  samples  used  in  this  series  of  experiments  consist 
of  zincblende,  chemical-vapor-deposition-grown  poly¬ 
crystalline  samples  of  ZnSe  and  CdTe  [5|.  The  ZnSe  sam¬ 
ple  was  3  mm  thick  and  the  CdTe  was  2  mm  thick. 


probe,  measured  outside  the  sample,  was  13°  and  the  peak 
input  irradiance  of  each  pump  was  lb  ~  34  MW  / cm2  and 
If  =  22  MW/cm2.  Clearly,  two  very  distinct  nonlineari¬ 
ties  are  evident  from  Fig.  2.  Near  zero  delay,  a  large  rap¬ 
idly  decaying  signal  is  seen,  while  at  longer  delays,  we 
observe  a  more  slowly  decaying  signal.  To  better  under¬ 
stand  the  two  nonlinear  regimes,  irradiance  dependence 
experiments  were  performed  at  different  delays.  Fig.  3 
shows  a  log-log  plot  of  the  DFWM  signal  versus  the  total 
input  irradiance,  (all  three  input  beams  were  varied  simul¬ 
taneously)  at  two  different  delay  times.  Fig.  3(a)  shows 
the  irradiance  dependence  at  zero  delay.  The  least-squares 
fit  gives  a  power  dependence  of  /3'2±0'2,  indicative  of  a 
third  order  nonlinearity  dominant  at  the  zero  delay  peak. 
Fig.  3(b)  shows  the  dependence  at  a  delay  of  240  ps,  with 
a  best  fit  giving  a  power  dependence  of  /5  0±0  2.  The  fifth 
order  dependence  of  the  DFWM  signal  on  the  input  beams 
can  be  explained  by  2PA  induced  carrier  refraction.  This 
mechanism  can  be  viewed  as  a  two  step  process.  First,  a 
modulated  carrier  density  is  generated  via  2PA;  this  is  an 
Im  {x(3)}  effect.  Second,  a  third  beam  diffracts  off  the 
carrier  modulation;  this  is  a  Re  {x(l1}  effect.  Hence,  the 
mechanism  is  a  sequential  x(3)  :x0)  process  that  appears 
as  an  overall  fifth  order  nonlinearity  [6].  Studies  in  CdTe 
at  1.06  /im,  where  this  material  exhibits  2PA,  reveal  the 
same  basic  behavior,  i.e,,  a  fast  third  order  signal  fol¬ 
lowed  by  a  slowly  decaying  fifth  order  signal. 


B.  Fast  Third  Order  Nonlinearity 

The  third  order  effect  observed  near  zero  delay,  as  can 
be  seen  from  Fig,  2,  decays  within  the  30  ps  pulsewidth, 
and  appears  unchanged  when  the  pump-probe  angle  is 
varied.  As  DFWM  is  sensitive  only  to  the  absolute  value 
of  x(3),  this  third  order  signal  may  arise  from  both  the  real 
and  imaginary  parts  of  the  susceptibility,  x{r  and  X/3\ 
where  xl3)  =  Xr  +  ‘X?'-  We  define  x(3)  in  terms  of  the 
nonlinear  polarization  in  e.g.s. -Gaussian  units  by;  P  = 
XtuE  +  x0)\E\2E/2.  For  all  beams  linearly  polarized 
parallel  to  each  other,  x<3)  =  6xim[ll],  as  we  are  dealing 
with  polycrystalline  (grain  size  =  1  /tm)  and  hence  effec¬ 
tively  isotropic  media.  The  source  of  x/3>  in  the  range 
Eg/ 2  <  hu  <  Eg,  is  2PA  (coefficient  02).  The  real  part 
of  x<3)  corresponds  to  nonlinear  refraction.  Usually,  this 
nonlinear  refraction  is  expressed  in  terms  of  n2,  where 
n2  =  n0  +  n2\E\2 /2.  n2  and  02  are  related  to  x<3>  by 
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n  2  =  —  Xr 
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(3) 
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III.  Results  and  Discussion 
A.  Identification  of  the  Nonlinear  Processes 

Using  pulses  at  0.532  /tm,  the  DFWM  signal  in  ZnSe 
was  monitored  as  a  function  of  input  energy  and  pulse 
delay  for  different  combinations  of  the  polarization  of  the 
three  input  beams.  Fig.  2  shows  a  plot  of  the  signal  versus 
the  delay  rb  of  Eh,  with  Eh  polarized  perpendicular  to  both 
Ef  and  Ep.  The  angle  9  between  the  forward  pump  and 
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where  e0  is  the  permittivity  of  free  space,  n0  is  the  linear 
refractive  index,  and  c  is  the  speed  of  light  in  vacuum. 
Here  n2  and  x<3)  are  expressed  in  esu,  and  all  other  param¬ 
eters  are  in  MKS  units.  We  can  estimate  the  magnitude 
of  n2  from  measurements  of  the  DFWM  reflectivity  tj  ep. 
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Fig.  2.  Phase-conjugate  signal  versus  backward  pump  delay  for  a  13° 
pump-probe  angle  and  with  the  backward  pump  polarized  perpendicular  to 
the  other  two  waves.  The  peak  pump  irradiances  are  If  =  22  MW  / cm2  and 
4  =  34  MW /cm2. 


IRRADIANCE  (MW/cm2) 

Fig.  3.  Log-log  plots  of  the  phase  conjugate  signal  versus  the  total  input 
irradianee  (If  +  Ih  +  Ip )  as  all  three  beams  were  varied  together,  for  (a) 
zero  delay  and  (b)  240  ps  delay.  The  solid  lines  are  best  fits  to  the  data 
giving  power  dependencies  of  73  1  ±0-2  and  /50±0'2,  respectively.  The  other 
experimental  parameters  are  as  in  Fig.  2. 


where  er  and  ep  are  the  energies  of  the  conjugate  and  probe 
pulses,  respectively.  In  the  degenerate  case  for  small  re¬ 
flectivity  and  in  the  “undepleted”  pumps  approximation 
the  third  order  susceptibility  is  given  by  [7] 


IxT  = 


n4c4e o  l 'c 

Mfh  l 


(3) 


where  /  is  the  peak  irradianee  for  each  beam,  n  is  the 
linear  refractive  index,  and  L  is  the  sample  thickness.  In 
(3),  the  ratio  of  irradiances  is  given  by 


(4) 


For  ZnSe  n  =  2.7  [9],  L  =  0.3  cm,  and  the  measured 
probe  and  calculated  conjugate  spot  radii  are  wp  —  1.41 
mm  (HWl  /e2M)  and  vv,  =  0.99  mm,  respectively.  tp/tc 
is  the  ratio  of  the  probe  to  conjugate  beam  pulsewidth 
which  for  a  third  order  effect  should  equal  V3,  giving  Ic/Ip 
=  0.069.  At  zero  delay  and  0.532  ftm  the  DFWM  energy 
reflectivity  of  ZnSe  was  measured  to  be  1.6%  for  If  =  46 
MW/cm2,  Ih  =  65  MW/cm2,  and  Ip  =  /,  /10.  The  2PA 
coefficient  for  ZnSe  at  this  wavelength  is  f}2  —  5.5 
cm/GW,  [10]  hence  at  these  irradianee  levels  losses  due 
to  2PA  are  less  than  6%,  so  that  the  undepleted  pumps 
approximation  is  valid.  This  gives  for  ZnSe: 

|X<3)|  =  2.6  x  10“19  m2/ V2  =  1.9  x  10"l2esu 

where  the  estimated  total  absolute  error  in  our  measure¬ 
ment  is  estimated  to  be  ±30%.  Substituting  /32  in  (2)  gives 
xf>  =  6.4  X  10“ 12  esu.  As  this  is  much  smaller  than  our 
measured  value,  this  indicates  that  there  must  also  be  a 
significant  real  component  to  x<3)-  Using 

lx(3>|  =  J(x?)2  +  (x?>f 

we  obtain  Ixr’I  =  1.8  x  10-12  esu.  Thus,  the  effect  of 
the  2PA  contribution  to  the  DFWM  signal  is  negligibly 
small,  given  our  estimated  errors  of  ±30%.  Expressing 
the  measured  Ixk'I  in  terms  of  n2  as  |n2|  —  4.2  X  10  11 
esu,  we  find  this  is  in  close  agreement  with  the  results  of 
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references  [1],  [8]  who  measured  n2  =  -4.4  x  10  1 1 
esu.  using  the  Z-scan  technique. 

In  order  to  further  investigate  the  nature  of  this  fast  third 
order  susceptibility,  we  also  performed  DFWM  experi¬ 
ments  in  polycrystalline  CdTe  at  1.064  p,m.  CdTe  ex¬ 
hibits  2PA  at  this  wavelength.  By  measuring  the  conju¬ 
gate  signal  in  CdTe  versus  the  delay  of  the  backward 
pump,  using  the  same  configuration  as  for  ZnSe  except 
using  1.064  /cm,  we  see  results  similar  to  those  shown  in 
Figs.  2  and  3  for  ZnSe.  At  zero  delay  we  find  that  the 
signal  varies  as  /31±0-3  while  at  long  delays  the  signal 
varies  as  /4  5±04.  For  CdTe  (3  =  22  cm/GW  [10]  and 
thus  beam  depletion  due  to  2PA  is  no  longer  negligible 
for  irradiance  levels  above  =50  MW /cm2.  The  some¬ 
what  reduced  slope  observed  at  long  delays  can  be  ac¬ 
counted  for  by  beam  depletion  at  higher  irradiances.  At 
an  irradiance  of  =  35  MW  / cm2  per  pump,  the  DFWM 
reflectivity  for  CdTe  is  14.8%  at  zero  delay.  This  yields 
a  third  order  susceptibility  of  |x<3>|  =  1.1  X  10  1 1  esu. 
The  2PA  coefficient  in  CdTe  at  1.06  g,m  has  been  mea¬ 
sured  as  j32  =  22  cm/GW.  [10]  Hence,  X/3)  =  5.6  x 
10“"  esu,  giving  jx/f’l  =  9.5  x  10“n  esu,  and  |n2|  = 
2.1  x  10“ 10  esu.  Once  again,  this  is  in  remarkably  close 
agreement  with  Z-scan  results,  which  gave  —2  X  10“ 10 
esu  [8], 

DFWM  experiments  were  also  performed  in  ZnSe  at 
1.064  (im  where  it  exhibits  three-photon  absorption.  At 
irradiance  levels  below  1  GW  /  cm2  per  pump  only  a  zero 
delay  peak  was  observed  that  varied  as  /2  8  *  0  3 ,  indicating 
the  dominance,  once  again,  of  a  third  order  nonlinearity. 
The  DFWM  signal  as  a  function  of  delay  of  the  backward 
pump  is  shown  in  Fig.  4,  while  in  Fig.  5  we  show  a  log- 
log  plot  of  the  irradiance  dependence  of  the  signal  at  zero 
delay.  As  the  phase-conjugate  reflectance  measured  at  this 
wavelength  Was  very  much  smaller  than  4%,  the  compar¬ 
ative  method  used  previously  for  determining  x <3>  could 
not  be  used  in  this  case.  Therefore,  x<3>  was  determined 
by  comparison  with  the  known  CS2  value  of  n2  =  1.3  X 
10“"  esu,  and  hence  Ixerrl  =  5.6  X  10“12  esu.  [12]  In 
order  to  make  such  a  comparative  measurement,  we  must 
account  for  differences  in  refractive  index,  surface  reflec¬ 
tance  and  lengths  of  the  two  samples  by  using  the  follow¬ 
ing  relationship  [13] 

I X  eff(ZnSe)  1 2 

lx$(CS2)|2 


nZnSe 


L  «cs. 


x  signal  ratio 


Lcsi  U  ~  files, 

£  ZnSe  (1  —  £)znSe 


(5) 


where  R  is  the  surface  reflectance  and  the  signal  ratio  is 
the  ratio  of  observed  conjugate  signal  energies  for  ZnSe 
and  CS2,  respectively.  Comparison  of  the  signals  from  the 
two  materials  under  identical  experimental  conditions 
gives  a  value  of  |x<3)|  =  1.2  ±  0.35  x  10“11  esu  for 
ZnSe  at  1.064  /xm.  As  there  is  no  2PA  in  this  wavelength 
region  for  ZnSe,  the  imaginary  part  of  x(3)  is  zero.  Thus, 
we  can  rewrite  our  result  as  |n2|  =  2.7  ±  0.9  x  10““ 


Backward  Pump  Delay  (ps) 


Fig.  4.  DFWM  signal  in  ZnSe  at  1.064  pm  versus  backward  pump  delay 
with  all  beams  parallel  and  both  pump  irradiances  —  600  MW/ cm2. 


Irradiance  (GW/cm2) 

Fig.  5.  Logarithmic  plot  of  the  DFWM  signal  in  ZnSe  at  1 .064  pm  versus 
the  total  incident  irradiance.  at  zero  delay.  The  solid  line  is  a  best  fit  to  the 
data,  giving  an  I2  8±03  dependence. 


esu.  This  also  agrees  within  errors  with  Z-scan  results, 
which  gave  n2  =  +1.7  x  10“u  esu  [8]. 


C.  Studies  of  the  Carrier  Nonlinearities  in  ZnSe  and 
CdTe 


Neglecting  diffusion  and  recombination  within  the 
picosecond  pulses,  [14]  the  growth  of  the  carrier  density 
N  generated  in  a  2PA  semiconductor  is  governed  by 


dN  ft,/2 
dt  2hw 


(6) 


where  (32  is  the  2PA  coefficient.  Here  we  have  ignored  the 
z  dependence  of  N  which  is  valid  for  our  experiments  in 
ZnSe  at  0.532  41m  where  the  irradiance  is  low  and,  there¬ 
fore,  the  2PA  is  small  ( <6%  loss).  In  the  DFWM  ge¬ 
ometry  / 2  =  (Ef  +  Eb  +  Epf .  We  concentrate  only  on 
those  terms  which  contain  EfE  *  (for  the  large  period  grat¬ 
ing)  and  neglect  terms  in  \Ep\2E*  since  \Ep\  «  \Ef\,  \Eb\. 
Then  the  component  of  /2  effective  in  producing  a  grating 
leading  to  the  DFWM  signal  is  /2ff  =  (n0ce0/2)2  (|£/|2  + 
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2\Eb\2)[2EfE*  +  c.c.].  Note  that  while  l\ ff  is  real,  it  has 
been  broken  down  into  the  sum  of  two  complex  conjugate 
terms.  The  first  gives  rise  to  the  conjugate  signal,  the  sec¬ 
ond  to  a  term  which  radiates  with  a  wave  vector  of  kp  — 
2k f,  and  hence  is  not  phasematched.  As  both  terms  are 
required  to  form  the  sinusoidal  grating,  we  must  retain 
both  of  them  to  calculate  the  solution  of  the  diffusion 
equation;  however,  the  effective  contribution  to  the  signal 
is  lift/ 2.  Thus,  the  effective  carrier  density  is  given  by 


The  resulting  nonlinear  polarization  source  term  PNL  is 
proportional  to  EhNclf.  In  the  case  where  Eh  is  not  tem¬ 
porally  overlapped  with  Ep  or  Ef,  /2ff  becomes  2(EfE*  + 
c.c. )\Ef  | 2  (effects  of  grating  decay  will  be  treated  in  the 
next  section).  Assuming  phase  distortions  on  Ep  to  be 
small,  E*  =  Epo  exp  {  +  ikp(z  cos  6  +  x  sin  0)}  and  Ef  = 
Ef0  exp  {—ikfz\ ,  so  that  can  be  written  as 


•  2  cos  [Kgx]  (  |£>(r')|3|£P(0|  dt'  (7) 

where  Kg  =  kf  —  kp  is  the  grating  wave  vector.  In  the 
limit  of  small  6,  Kg  —  2x6 /X,  where  X  is  the  free-space 
wavelength.  Since  Ec  oc  />NL,  the  total  irradiance  depen¬ 
dence  of  the  DFWM  signal  is 

h  <*  hi Pl}-  (8) 

In  order  to  verify  this,  we  measured  the  DFWM  signal 
in  CdTe  at  1.064  p m  as  a  function  of  /,,  only  with  the 
backward  pump  delayed  by  240  ps.  The  signal  was  plot¬ 
ted  logarithmically  against  and  the  slope  of  the  resulting 
straight  line  was  found  to  be  2.8  ±  0.3  in  accordance  with 
(8). 

A  similar  expression  can  be  obtained  when  the  carriers 
are  generated  via  a  three-photon  absorption  process.  Fol¬ 
lowing  an  analogous  procedure,  the  end  result  is  a  seventh 
order  dependence  of  the  DFWM  signal  on  the  input  beams 
as  three-photon  absorption  is  an  Im  x<5)  process.  In  our 
experiments  in  ZnSe  at  1.064  /tm,  we  find  that  at  higher 
input  irradiances,  on  the  order  of  1 .2  GW / cm2  per  strong 
pump,  a  free-carrier  tail  became  evident.  Fig.  6  shows  a 
log-log  plot  of  the  irradiance  dependence  for  rb  =  240 
ps.  A  best  fit  gives  a  dependence  of  /6'9±n2,  agreeing 
closely  with  the  power  of  7  expected  for  3  PA  induced 
carrier  refraction.  Under  our  experimental  conditions  with 
the  backward  pump  delayed,  the  DFWM  signal  depen¬ 
dence  is 

Ic  oc  lblpl).  (9) 

Fig.  7  shows  a  similar  plot  to  Fig.  6,  except  that  in  this 
case  only  If  is  varied  resulting  in  a  power  dependence  of 
Is  1  ±0  3,  close  to  the  power  of  5  predicted  in  (9). 


Pump  Irradiance  (GW/cm2) 

Fig.  6.  Logarithmic  plot  of  the  DFWM  signal  resulting  from  the  three- 
photon  absorption  generated  carrier  grating  in  ZnSe  at  1 .064  /zm  versus  the 
total  input  irradiance.  The  measurement  was  performed  with  the  backward 
pump  delayed  by  240  ps.  The  solid  line  is  a  best  fit  to  the  data,  indicating 
an  /69±02  dependence. 
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Fig.  7.  As  Fig.  6,  but  with  only  the  forward  pump  varied  and  the  backward 
pump  and  probe  held  at  constant  irradiance.  The  forward  pump  energy  was 
not  calibrated  for  this  experiment.  The  solid  line  is  the  best  fit,  giving  an 
1'j  ±0 3  dependence  for  the  DFWM  signal  in  this  case. 


D.  Temporal  Evolution  of  Diffraction  Efficiency  in  ZnSe 

In  order  to  determine  the  refractive  index  change  per 
excited  carrier  pair  per  unit  volume,  on  (units  of  cm3),  we 
need  to  measure  the  diffraction  efficiency  of  the  carrier 
grating  in  the  absence  of  the  third  order  signal.  This  ne¬ 
cessitates  measurement  of  the  DFWM  signal  at  a  delay 
long  compared  to  the  pulsewidth.  a„  is  thus  obtained  by 
fitting  the  theoretical  diffraction  efficiency  to  the  mea¬ 
sured  value.  However,  to  calculate  this  diffraction  effi¬ 
ciency  we  require  a  knowledge  of  the  temporal  depen¬ 
dence  of  the  density  of  2PA-generated  carriers  N(x,  t). 
This  is  governed  by  [15] 

dNcff(x,  t )  ■  Ncff(x,  t)  _  dzNett(x,  t)  =  g/2ff(x,  t) 

dt  t  dx2  2hw 

GO) 
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where  r  is  the  carrier  decay  time  and  D  is  the  carrier  dif¬ 
fusion  coefficient.  We  use  a  constant  carrier  lifetime  as 
the  carrier  densities  in  our  experiments  are  several  orders 
of  magnitude  too  low  to  have  significant  bimolecular  or 
Auger  recombination.  This  is  verified  by  the  Is  power  de¬ 
pendence  in  our  experiments.  Since  the  beam  radius  is 
much  larger  than  the  grating  period  we  assume  the  diffu¬ 
sion  occurs  only  in  the  direction  x  between  grating  max¬ 
ima  and  minima.  Further,  we  have  assumed  low  excita¬ 
tion  such  that  the  spatial  amplitude  of  I2  does  not  vary  in 
the  propagation  direction  z,  i.e.,  j3 IL  «  1  which  was 
true  for  the  irradiance  levels  used.  The  solution  to  the  dif¬ 
fusion  equation  is  then 

ft 

Neft(x,  t)  =  — -  exp  (-t/Tg) 


exp  (t'/rg)ll ff(x,  y,  t')  dt’ 


(11) 


where  jg  is  the  total  decay  time  of  the  grating,  given  by 
Tg  =  tdtr/(td  +  tr)-  Here  td  is  the  grating  diffusion 
lifetime  and  tr  is  the  carrier  recombination  lifetime  [14], 
We  have  again  assumed  that  Eb  is  time  delayed.  Note  that 
the  product  EfE*  in  l\ ff  has  a  term  2\Ef\\Ep\  exp  (iKgx), 
where  Kg  —  2ir /A,  and  A  =  X / [2 n  sin  (0/2)]  is  the  grat¬ 
ing  spacing.  Hence  l/rs  =  (4ir2D)/A2  +  1/r  [16].  For 
the  sample  thickness  and  grating  spacings  used,  we  were 
always  in  the  thick  grating  limit.  The  diffraction  effi¬ 
ciency  for  weak  excitation  of  a  thick  lossless  phase  grat¬ 
ing  was  derived  by  Kogelnik  [4],  [17]  to  be 


7)  -  sin2 


/  AnirL  \ 
\X  cos  8/ 


'  AnirL  \' 
A  cos  6J 


(12) 


where  An  =  <r„Afcfr.  The  energy  diffracted  into  the  con¬ 
jugate  beam  et.  is  then  given  by  temporally  and  spatially 
integrating  the  product  of  the  diffraction  efficiency  rj  and 
the  beam  irradiance  incident  on  the  grating  giving 

foo 

n(x',  y',  t')Ib(x',  y',  t'  -  Tb)  dx'  dy'  dt’ 

—  00 


(13) 

where  Ib(x,  y,  t)  =  Ib0  exp  [-(.r2  +  y2)/w2  -  (f2/^)2]. 
Here  rb  accounts  for  the  time  delay  of  the  backward  pump 
with  respect  to  the  interfering  beams.  Note  for  t/  =  1 ,  (13) 
gives  the  energy  in  the  backward  pump  eh.  Defining  y,  as 
the  energy  diffraction  efficiency,  we  find 


V( 


ec/ Eft  —  £(1  +  £)' 


2  V2ff  2<r  n  //  "Jxtp  L 7 
80  (htdfX2  cos2  d 


exp  (c/4t2) 


Joo 

exp  (-(*'  “  r'>)2/^)  exp  (-2 t'/Tg) 


•  {1  +  erf  [^2/tp(t'  -  t2p/ 4rg)]}2  dt'  (14) 

where  t,  is  the  ratio  of  lp  to  /.  Note  that  this  derivation  is 
valid  only  when  the  beam  widths  are  large  compared  to 
the  grating  period,  when  the  length  of  interaction  of  the 
interfering  beams  is  much  larger  than  the  sample  thick- 
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V) 


Carrier  Density  (xIO15  cm'3) 

Fig.  8.  Plot  of  measured  diffusion  coefficient  versus  peak  carrier  density. 


ness,  and  when  absorption  of  the  interacting  beams  is 
small.  Given  ft  =  5.5  cm/GW  [10]  at  532  nm  for  ZnSe, 
a  numerical  fit  of  (14)  yields  a„  and  rg. 

The  experimentally  measured  value  for  the  diffraction 
efficiency  for  ZnSe  at  0.532  pm  at  r  =  240  ps  is  tj  =  6.8 
x  10  5  for  a  pump-probe  angle  6  of  13°.  The  grating 
decay  rate  rg  was  determined  from  the  best  fit  exponential 
decay  to  the  data  at  long  delays  for  several  input  irradi- 
ances.  Using  tp  =  18  ps  (HW1  /eM  in  irradiance),  and  Tg 
=  247  ps,  and  numerically  integrating  (14),  we  obtain 
!<j„  =5.1  +2.5  x  10  1  cm1.  This  is  considerably  larger 
than  the  result  of  another  recent  measurement  of  o„  that 
obtained  a„  =  -0.8  +  0.2  x  10~21  [18]  and  probably 
indicates  the  difficulty  of  obtaining  accurate  absolute  val¬ 
ues  of  nonlinearities  from  multiple  beam  diffraction  effi¬ 
ciency  measurements.  Note  that  the  measurement  in  [18] 
showed  the  sign  of  a„  to  be  negative  [3],  The  uncertainties 
stated  are  relative  fitting  errors.  The  absolute  error  bars 
on  tr„  are  estimated  to  be  =50%. 

We  were  also  able  to  extract  the  density  dependence  of 
the  carrier  diffusion  coefficient  by  measuring  the  angular 
dependence  of  the  grating  decay  time  for  a  number  of  dif¬ 
ferent  pump  irradiances  and  therefore  different  carrier 
densities  [19].  This  dependence  is  plotted  in  Fig.  8.  We 
find  that  for  irradiances  above  40  MW /cm2,  correspond¬ 
ing  to  carrier  densities  >5  x  1014  cm-3,  D  maintains  a 
constant  value  of  4.5  +  0.5  cm2/s,  while  below  5  x  1014 
cm'/  the  diffusion  coefficient  decreases  rapidly.  For  these 
densities  we  do  not  expect  to  have  any  density  depen¬ 
dence  of  either  the  mobility  or  the  recombination  lifetime. 
We  therefore  propose  that  the  small  diffusion  coefficient 
at  low  density  may  be  due  to  trapping  of  the  2PA  excited 
carriers.  As  the  carrier  density  is  increased,  these  traps 
become  filled  and  the  diffusion  is  dominated  by  the  am- 
bipolar  mobility.  This  high  density  result  is  comparable 
to  the  result  of  2.5  +  0.5  cm2/s  obtained  by  Jarasiunas 
and  Gerritsen,  but  their  carrier  density  is  not  given  [3], 
Using  data  from  Hall  mobility  measurements  in  ZnSe, 
[20]  we  calculate  an  ambipolar  mobility  of  26  cm2  V  1 
s-,  corresponding  to  a  diffusion  coefficient  of  0.7  cm2/s 
at  300°K.  This  is  in  good  agreement  with  our  low  density 
measurements,  since  the  Hall  mobility  measurements  are 
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performed  with  low  dopant  densities.  It  should  be  noted 
that  over  the  entire  range  of  carrier  densities  measured, 
the  carrier  nonlinearity  is  fifth  order.  This  rules  out  the 
possibility  that  the  traps  are  filled  by  linear  absorption. 
These  measurements  also  gave  an  estimate  for  the  recom¬ 
bination  lifetime  of  tr  =  725  ±  275  ps. 

IV.  Conclusion 

We  have  presented  a  series  of  picosecond  transient 
DFWM  measurements  performed  in  two  direct-gap  II- VI 
semiconductors.  For  ZnSe,  large  third  order  nonlinear 
mechanisms  were  observed  both  at  0.532  and  1.064  fim. 
The  origin  of  these  fast  third  order  nonlinearities  is  be¬ 
lieved  to  be  due  to  both  two-photon  absorption  and  the 
bound  electronic  Kerr  effect  [8],  The  magnitude  of  the 
nonlinear  susceptibilities  found  were  |x<3,|  —  1.9  ±  0.57 
x  10-11  esu  and  1.2  +  0.35  x  10“"  esu  at  0.532  /cm 
and  1.064  /on,  respectively.  Similarly,  for  CdTe  a  third 
order  effect  yielded  a  |x'3,|  =  7.2  +2.1  X  10“ 10  esu  at 
1.064  fi m.  Combined  with  independent  measurements 
[10]  of  the  two-photon  absorption  coefficients  in  these 
materials,  this  indicates  that  the  measured  susceptibility 
in  ZnSe  at  532  nm  and  CdTe  at  1.06  fit n  is  dominated  by 
nonlinear  refraction.  As  there  is  no  2PA  in  ZnSe  at  1 .06 
fim,  we  conclude  that  this  susceptibility  is  entirely  refrac¬ 
tive.  Other  data  [8]  indicates  that  in  ZnSe  n2  is  negative 
at  0.532  fim  and  positive  at  1 .064  /on,  while  in  CdTe  at 
1.064  /on,  n2  is  negative.  Higher  order  nonlinearities 
(higher  than  third  order)  were  also  observed  in  both  semi¬ 
conductors.  These  nonlinearities  are  attributed  to  nonlin¬ 
ear  absorption  induced  carrier  refraction. 

When  carriers  are  excited  via  the  absorption  of  two 
photons,  we  have  seen  a  fifth  order  dependence  of  the 
conjugate  field  on  the  input  fields  as  expected.  In  the  same 
manner,  for  carrier  excitation  achieved  via  the  absorption 
of  three  photons,  an  effective  seventh  order  nonlinearity 
was  found. 

For  ZnSe  at  0.532  /ira  and  for  low  excitation,  a  value 
for  the  index  of  refraction  per  carrier  pair  generated  of  a„ 
=  5.1  +  2.5  X  10  21  cm3  was  estimated.  By  monitoring 
the  grating  decay  times  for  different  grating  spacings,  the 
density  dependence  of  the  diffusion  coefficient  was  ex¬ 
tracted. 

Acknowledgment 

We  wish  to  thank  B.  S.  Wherrett  for  his  critical  review 
of  this  paper,  G.  M,  Moharam  for  useful  discussions  on 
the  diffraction  efficiency  calculations,  and  E.  Miesak  for 
help  in  taking  data. 

References 

fl]  M.  Sheik-Bahae,  A.  A.  Said,  T.  H.  Wei.  D.  J.  Hagan,  and  E.  W. 
Van  Stryland,  "Sensitive  measurement  of  optical  nonlinearities  using 
a  single  beam,”  IEEE  J.  Quantum  Electron.,  vol.  26,  pp.  760  769, 
1990. 


[2]  R  K.  Jain  and  M  B.  Klein,  “Degenerate  four-wave  mixing  in  semi¬ 
conductors,’'  in  Optical  Phase  Conjugation ,  R.  A.  Fisher.  Ed.  New 
York:  Academic,  1983.  pp.  307-415. 

[3]  K.  Jarasiunas  and  H.  J.  Gerritsen,  “Ambipolar  diffusion  measure¬ 
ments  in  semiconductors  using  nonlinear  transient  gratings."  Appl. 
Phvs.  Lett vol.  33,  pp.  190-193,  1978. 

{4]  H.  Kogelnik,  “Coupled  wave  theory  for  thick  hologram  gratings," 
Bel!  Syst.  Tech.  J.,  vol.  48,  pp.  2909-2947,  1969. 

[5]  Sample  obtained  from  II— VI.  Inc.,  Saxonburg,  PA. 

|6]  E.  W.  Van  Stryland,  Y.  Y.  Wu.  D.  J.  Hagan,  M.  J.  Soileau,  and  K. 
Mansour,  “Optical  limiting  with  semiconductors."  J.  Opi.  Soc. 
Amer.  B.,  vol.  5,  pp.  1980-1989,  1988. 

[7]  J.  J.  Wynne,  "Optical  third-order  mixing  in  GaAs,  Ge,  Si,  and  InAs," 
Phys.  Rev.,  vol.  178,  pp.  1295-1303,  1969. 

]8]  M.  Sheik-Bahae,  D.  J.  Hagan,  and  E.  W,  Van  Stryland,  “Dispersion 
and  band-gap  scaling  of  the  electronic  Kerr  effect  in  solids  associated 
with  two-photon  absorption,"  Phys.  Rev.  Lett.,  vol.  65,  pp.  96-99, 
1989. 

|9]  Landolt-Bornstein,  Numerical  Data  and  Functional  Relationships  in 
Science  and  Technology ,  Volume  17,  Semiconductors,  K.-H.  Hel- 
lewege,  Ed.  New  York:  Springer- Verlag.  1982. 

(10J  E.  W.  Van  Stryland,  H.  Vanherzeele,  M.  A.  Woodall.  M.  J.  Soileau, 
A.  L.  Smirl,  S.  Guha,  T.  F.  Boggess.  “Two-photon  absorption,  non¬ 
linear  refraction  and  optical  limiting  in  semiconductors."  Opt.  Eng., 
vol.  24,  pp.  613-623,  1985. 

Ill]  P.  N.  Butcher.  “Nonlinear  optical  phenomena."  Bulletin  200,  Ohio 
State  Univ.,  Columbus.  OH,  1965. 

[12]  M.  J.  Soileau,  W.  E.  Williams,  and  E.  W.  Van  Stryland,  “Optical 
power  limiter  with  picosecond  response  time.”  IEEE  J.  Quantum 
Electron.,  vol.  QE-19,  p.  731.  1983. 

[13]  R.  Adair,  L.  L.  Chase,  and  S.  A.  Payne,  “Nonlinear  refractive  index 
of  optica!  crystals,”  Phys.  Rev.  B,  vol.  39.  pp.  3377-3350,  1987. 

[14]  J.  H.  Bechtel  and  W.  L.  Smith,  “Two-photon  absorption  in  semicon¬ 
ductors  with  picosecond  laser  pulses,"  Phys.  Rev.  B.  vol.  13,  pp. 
3515-3522,  1976. 

[15]  W.  Van  Roosbroeck,  “The  transport  of  added  current  carriers  in  a 
homogeneous  semiconductor,”  Phvs.  Rev.,  vol.  91.  pp.  282-289, 
1953. 

[16]  J.  R.  Salcedo,  A.  E.  Siegman,  D.  D.  Dlott,  and  M.  D.  Fayer,  “Dy¬ 
namics  of  energy  transport  in  molecular  crystals:  the  picosecond  tran¬ 
sient  grating  method,"  Phys.  Rev.  Lett.,  vol.  41,  p.  131,  1978. 

[17]  H  J.  Eichler,  P.  Gunter,  and  D.  W.  Pohl.  Laser-Induced  Dynamic 
Gratings,  Springer  Series  in  Optical  Sciences,  J.  M.  Enoch,  D.  L. 
MacAdam,  A.  L.  Schawlow,  K.  Shimoda,  and  T.  Tamir.  Eds.  New 
York:  Springer-Verlag,  1988. 

[18]  A.  A.  Said,  M.  Sheik-Bahae,  D.  J.  Hagan,  E.  J.  Canto-Said,  Y.  Y. 
Wu,  J.  Young,  T.  H.  Wei.  and  E.  W.  Van  Stryland,  “Nonlinearities 
in  semiconductors  for  optical  limiting.”  in  Proc.  SPIE,  Electro- 
Optical  Materials  for  Switches.  Coatings.  Sensor  Optics,  and  Detec¬ 
tors,  Orlando,  1990,  conf.  1370. 

[19]  D.  J.  Hagan,  H.  A.  MacKenzie.  H.  A.  Al  Attar,  and  W.  J.  Firth, 
“Carrier  diffusion  measurements  in  InSb  by  the  angular  dependence 
of  degenerate  four-wave  mixing.”  Opt.  Lett.,  vol.  10,  p.  187,  1985. 

[20]  B,  Ray.  1 1 -VI  compounds  .  Edinburgh:  Pergamon.  1969,  ch.6. 


E.  J.  Canto-Said,  photograph  and  biography  not  available  at  the  time  of 
publication. 


D.  J.  Hagan  (M’87),  for  a  photograph  and  biography,  see  p.  1309  of  the 
June  1991  issue  of  this  Journal. 


J.  Young,  photograph  and  biography  not  available  at  the  time  of  publica¬ 
tion. 


Eric  \V,  Van  Stryland  (M’84-SM’90).  for  a  photograph  and  biography, 
see  p.  1309  of  the  June  1991  issue  of  this  Journal. 


